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Abstract 

It is demonstrated that in low multiplicity sample, the increase 
of the fluctuation of event-factorial-moments with the diminishing of 
phase space scale, called "erraticity" , are dominated by the statistical 
fluctuations. The erraticity behavior observed at NA27 experiment 
can be readily reproduced by pure statistical fluctuations. Applying 
erraticity analysis to a high multiplicity sample is recommended and 
the method is improved at very high multiplicity well. 
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Since the finding of unexpectedly large local fluctuations in a high multiplicity event 
recorded by the JACEE collaboration |l|], the investigation of non-linear phenomena in high 
energy collisions has attracted much attention ||. The anomalous scaling of factorial mo- 
ments, defined as 

p 1 ^4 (n m (n m -!)■■■ (n m - q + 1)) 

at diminishing phase space scale or increasing division number M of phase space ||: 

F q oc M _< \ (2) 

called intermittency (or fractal) has been proposed for this purpose in multiparticle system. 
The average (■ • •) in Eqn.(l) is over the whole event sample and n m is the number of particle 
falling in the mth bin. That kind of anomalous scaling has been observed successfully in 
various experiments @]||. 

A recent new development further along that direction is the event-by-event analy- 
sis 00- An important step in this kind of analysis was made by Cao and Hwa [Q, 
who first pointed out the difference between a dynamic system in which time sequence can 
be traced and the multiparticle system where only (phase) space patterns can be obtained. 
They proposed to measure the pattern by the event factorial moments 
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as opposed to sample-factorial-moments defined in Eqn.(l) averaged over all events. Its 
fluctuations from event to event can be quantified by its normalized moments as: 

<? M =<*?>, % = F^/(F^), (4) 

and by dC Pjq /dp at p = 1: 

Z q = (%\n%) (5) 

If there is a power law behavior of the fluctuation as division number goes to infinity, or as 
resolution 5 = A/M goes to very small, i.e., 

C M {M) oc M^ (p) , (6) 



the phenomenon is referred to as erraticity 0. The derivative of exponent ip q (p) at p — 1 



describes the width of the fluctuation and so is called as entropy index. In the following, we 
will call C Ptq or £ g in Eqn.(4) and (5) as erraticity-moments. 

It is well known that the obstacle of event-by-event analysis is the influence of statistical 
fluctuations caused by insufficient number of particles. The big advantage of sample factorial 
moments in Eqn.(l) is that it can eliminate this kind of statistical fluctuations. It has been 
proved || that if the statistical fluctuations of particles falling in a bin is Possonian like, then 
the sample-factorial-moments equal to the corresponding dynamic probability moments: 

M 

F q = C q = M^J2(pl)- (8) 
1=1 

Again, the average is over the whole sample. The pi is the probability of particle falling in 
the ith bin in a certain event. However, we can not follow the same procedure to get the 
similar equation for event factorial moments of Eqn.(3). Since the number of particles in 
an event is not large enough and so is the number of bin, event factorial moments can not 
completely eliminate statistical fluctuations and present the dynamic probability moments 
associated with it. How large of the statistical fluctuations in erraticity analysis is and how 
it depends on the number of multiplicity have not been seriously estimated yet. We will 
answer these questions quantitatively in the letter. 

To be direct and obvious, we firstly use an unique flat probability distribution in whole 
studying interval and whole sample. It means that the probabilities in all bins are equal and 
are the same for different events. For simplicity, we use only the fixed number of multiplicity. 
In this case the denominator in the definition of factorial moment, eqn.(l) and (3), becomes 
simply N(N — 1) ||. We first take N = 9, which is about the average multiplicity at ISR 
energies. The distribution of particle in the whole studying phase region in an event can be 
readily located by Bernouli distribution: 

AT! M 

B(n 1 ,...,n M \p 1 ,...,p M ) = — 7Pi 1 ---PM t i z2 n m = N. (9) 

By this way, we simulate a sample with 1000 events. The results of the second order sample- 
factorial- moments F2, erraticity-moments C Pt 2 and £2 on the division number M of the 
phase-space region are shown in Fig. 1(a). The second order sample-factorial- moment is a 
constant with the increasing of division number. This is what we expect. Since no dynamic 
is input, it becomes a constant after eliminating the statistical fluctuations. While the 
increase of erraticity-moments C Pt 2 and £2 with division number is measurably large. These 
contributions come from pure statistical fluctuations of event factorial moments due to the 
insufficient number of particle in an event since there is no dynamic fluctuation from event 
to event in the case. This results can fully recover what has observed in NA27 data [HJ, cf. 



the open circles in the second figure of Fig. 1(a). This means that in low multiplicity sample, 
the statistical fluctuations of event factorial moments dominate the erraticity behavior of 



multiparticle system. Event factorial moments is not a good representation of event dynamic 
at low multiplicity events. 

However, erraticity analysis proposed a very important way to study the event-by-event 
fluctuations though we are still not clear whether there are such fluctuations and if there are 
what mechanism causes them. We have demonstrated in our former paper [jn| that if and 
only if different events have different dynamic fluctuation strengths, the erraticity moments 
will keep increasing with the increasing of division number and so has nonzero entropy index. 

As is well known, the statistical fluctuations will become negligible if the multiplicity 
of an event is large enough. At how high a multiplicity the event factorial moments can 
measure the dynamic fluctuations of a finite particle system is a very meaningful question. 
In the left of the paper, we will focus our discussion on answering the question. 

Now we switch the fixed multiplicity iV to 20 and 300 in the above mentioned simulation. 
The corresponding second order sample-factorial-moments F2, erraticity- moments C p ^ and 
S 2 versus the division number M of the phase-space region are shown in Fig. 1(6) and (c) 
respectively. The second order sample-factorial-moment keep to be a constant as we know. 
The erraticity moments become flater and flater with the increase of multiplicity. It means 
that pure statistic fluctuations of event factorial moments are greatly depressed by the 
increase of multiplicity. 

From Fig.(l), we can see that, when multiplicity is larger than 300, event factorial mo- 
ments can be approximately used to describe the event spatial pattern associated with it and 
its moments — erraticity moments — can represent the erraticity behavior of the system 
safely. 

In order to confirm this upper limit of multiplicity, we do following parallel analysis for 
a system with dynamic fluctuation from event to event. As we know ||l2"f , random cascading 
model, or a-model is the simplest model which can be used to generate a sample with nonzero 
entropy index. We will use it for our quantitative discussion below. In the random cascading 
a-model, the M division of a phase space region A is made in steps. At the first step, it is 
divided into two equal parts; at the second step, each part in the first step is further divided 
into two equal parts, and so on. The steps are repeated until M = AY/5y = 2 U . How 
particles are distributed from step-to-step between the two parts of a given phase space cell 
is defined by the independent random variable u) V j v , where j u is the position of the sub-cell 
(1 < jv < 2") and v is the number of steps. It is given by flTTf : 

I 1 
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where, r is a random number distributed uniformly in the interval [—1,1]. a is a positive 
number less than unity, which determines the region of the random variable us and describes 
the strength of dynamical fluctuations in the model. If it change from event to event, there 
will be different dynamic fluctuation strength in different events. Here, let it has a Gaussian 
distribution. The mean and variance of the Gaussian are both chosen as 0.22. After v steps, 



the probability in the rath window (m = 1, . . . , M) is p m = wij 1 U2j 2 ■ ■ ■ oo V j v . 

By the model, we generate 1000 events. The intermittency analysis, or the logarithm 
of second order sample probability moment lnC2, and erraticity moments lnC Pi 2 and E2 
as function of InM are shown in Fig.2(a). Now the sample probability moment MC2 has 
a power law behavior as dynamic fluctuations have been input. The erraticity moments 
also show a power law behavior at large division number region. It represents the dynamic 
fluctuation from event to event. The corresponding entropy index obtained from a linear fit 
to the last 3 points of £ 2 is /i 2 = 0.0161. 

Finite number of particle can also be added to the above pure dynamic fluctuation model 
by Bernouli distribution of Eqn.(9). Again, we put N = 9 first. The corresponding factorial 
moment and erraticity moments are given in Fig.2(c). The value of erraticity moments now 
are much larger than those obtained from the original pure dynamic fluctuations in Fig.2(a). 
The entropy index, ^2 = 0.422, also turns out to be more than one magnitude bigger. This 
results confirm us again that the erraticity behavior is dominated by statistical fluctuations 
in low multiplicity sample if we use event factorial moments to characterize it. Though there 
is dynamic fluctuation from event to event, it will be merged to large statistical fluctuations 
in the case. 

Secondly, we let iV = 300. The corresponding factorial and erraticity moments are shown 
in Fig. 2 (b). The erraticity moments now approach to its original dynamic fluctuation values 
in Fig. 2(a) and entropy index is fi 2 = 0.0168 close to its real value 0.0161. So we get the same 
conclusion as pure statistic fluctuation case. After multiplicity is larger than 300, erraticity 
behavior of the system can be estimated by the fluctuation of event factorial moments. 

To show quantitatively the influences of statistic fluctuations on erraticity behvior at 
different multiplicity cases, we simulate various number N = 5, . . . , 1000 of particles in 
an event for both flat probability distribution and above-described a model cases. The 
corresponding entropy indices are given in Fig. 3 as full circles (flat distribution) and full 
triangles (a model) respectively. We can see that both of them decrease with multiplicity. 
For flat probability distribution, entropy index of statistical fluctuation is depressed more 
than three orders of magnitude when multiplicity N increases from a few to 300. After 
multiplicity N > 300, the entropy index is unmeasurably small. Meanwhile, the entropy 
index of a-model sample approaches to its real dynamic value /i 2 = 0.0161, cf. the solid line 
in Fig. 2, after N > 300. The multiplicity of current and nearly coming heavy-ion collision 
is about this number or higher. The erraticity analysis given by event factorial moments is 
applicable for heavy-ion collisions, when the average multiplicity is higher than 300, where 
we are free from the influence of the statistical fluctuations. 

In fact, if multiplicity is higher than a thousand, which has been recorded in NA49 ex- 
periments and will be the case for the future heavy-ion collision experiments, the factorial 
moments analysis of a single event is unnecessary anymore as n m (n m — 1) • • • does not make 
much difference from n m ■ n m • • • in most of the phase-space bins which provides main con- 



tribution to the anomalous scaling of moments. In these cases the probability distribution 
in an event can be approximately presented by: 



n 



M 



Pm 



rn 



1 = H Pm- 



(11) 



m=l 



The erraticity-moments of event-probability-moments can be consequently defined by: 



By this definition, we repeat the analysis for both the flat probability distribution and the 
dynamic-fluctuation distribution cases. It is a little bit smaller than the corresponding event- 
factorial-moments analysis at flat probability distribution case, cf. the open and full circles 
in Fig .3, and so it depresses the influence of statistical fluctuations more. The dynamic 
fluctuations of event to event in the a model case can still be abstracted out as done by the 
event-factorial-moments description, cf. the open and full triangles in Fig. 3. 

From the simple discussion above, we can make the following conclusions: If we use event 
factorial moments to measure spatial pattern, in very low multiplicity sample, such as the 
sample of ISR energies, statistical fluctuations caused by insufficient number of particle in an 
event will control the erraticity behavior of the system. Therefore, the physical conclusions 
from the experimental data on these kind of sample can not be treated seriously. However, if 
the multiplicity of studying sample is larger than 300, the influence of statistical fluctuations 
on erraticity behavior is negligible. Therefore, the erraticity behavior, if any, could be well 
observed in current and future heavy-ion collisions. Further more, if the multiplicity of an 
event is larger than a thousand, the probability moments defined by Eqn.(ll)(12) can present 
erraticity behavior of the system as well as the event-factorial-moments. 
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Figure Captions 



Fig.l (a) The dependence of the logarithm of the second order sample-factorial- 
moments F 2 , erraticity-moments C Pi2 and £ 2 on that of the phase-space 
division number M for a flat probability distribution with particle number 
equal to 9 (a), 20 (6) and 300 (c) respectively. The dashed lines are linear 
fit. Open points are experimental results of NA27. The solid points are MC 
results. The solid lines are for guiding the eye. 

Fig. 2 (a) The dependence of the logarithm of the second order probability- moments 
C 2 , erraticity-moments C Pi2 and £ 2 on that of the phase-space division 
number M for the a model with Gaussian-distributed a. (b) The same as (a) 
but for the sample-factorial-moments F 2 and the corrssponding erraticity- 
moments C Pt 2 and S 2 with 300 particles, (c) The same as (b) but with 9 
particles. The dashed lines are linear fit. The solid lines are for guiding the 
eye. 

Fig. 3 The dependence on number of particle of the entropy indices /i 2 for Gaussian 
a-model calculated from event-factorial-moments (full triangles) and from 
probability- moments (open triangles). The same for flat distribution (full 
and open circles). The solid line is the dynamical result without statistical 
fluctuation. The dashed lines are for guiding the eye. 
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Fig.l (a) The dependence of the logarithm of the second order sample-factorial- 
moments F 2 , erraticity- moments C Pi2 and S 2 on that of the phase-space 
division number M for a flat probability distribution with particle number 
equal to 9 (a), 20 (b) and 300 (c) respectively. The dashed lines are linear 
fit. The solid lines are for guiding the eye. 
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Fig. 2 (a) The dependence of the logarithm of the second order probability- 
moments C*2, erraticity- moments C, p ^ and S 2 on that of the 
phase-space division number M for the a model with Gaussian- 
distributed a. (b) The same as (a) but for the sample-factorial- 
moments i<2 and the corrssponding erraticity-moments C Pi 2 and 
S 2 with 300 particles, (c) The same as (b) but with 9 particles. 
The dashed lines are linear fit. The solid lines are for guiding 
the eye. 
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Fig. 3 The dependence on number of particle of the entropy indices [12 
for Gaussian a-model calculated from event-factorial-moments 
(full triangles) and from probability- moments (open triangles). 
The same for flat distribution (full and open circles). The solid 
line is the dynamical result without statistical fluctuation. The 
dashed lines are for guiding the eye. 



